Quantum fluctuations in one-dimensional arrays of condensates 
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The effects of quantum and thermal fluctuations upon the fringe structure predicted to be observable 
in the momentum distribution of coupled Bose-Einstein condensates are studied by the effective- 
potential method. For a double- well trap, the coherence factor recently introduced by Pitaevskii and 
Stringari [Phys. Rev. Lett. 87, 180402 (2001)] is calculated using the effective potential approach 
and is found in good agreement with their result. The calculations are extended to the case of a one- 
dimensional array of condensates, showing that quantum effects are essentially described through 
a simple renormalization of the energy scale in the classical analytical expression for the fringe 
structure. The consequences for the experimental observability are discussed. 



Italy, 



Recently, extensive theoretical work was devoted to the 
problem of the coherence properties arising from the in- 
teraction among two or more trapped condensates |Q-j§] . 
These effects are experimentally accessible and can 
be considered of the same nature as in Josephson junc- 
tions ||. As a matter of fact, the condensates can 
be characterized by their phases 0i and atom numbers 
hi, which are conjugate canonical variables satisfying 



[0i,'< 



iSi-j. The Hamiltonian looks indeed like the 



Josephson one 



namely 



n = 



E r 



E 

i=\ 



E 



(ij) 



<Pj): 



(1) 



for an array of N s condensates, where the second sum- 
mation runs over the nearest-neighbor condensates. For 
two interacting condensates (i = 1, 2) this reduces to 



Tin 



E,. 



E 3 cos 



(2) 



01 ~ 02, and the unessential 
-7T-2) 2 is omitted. The inter- 



where h = (hi — fi2)/2, = 
conserved term (E^/8) (fr- 
action parameters E c and E 1 depend on the chemical 
potential and the wave-function overlap of neighboring 
condensates, respectively ||. For the above Hamiltoni- 
ans it is convenient to define a quantum coupling pa- 
rameter g as the ratio between the quantum energy scale 
hujo = \J E c E 3 of the quasiharmonic excitations and the 
overall energy scale of the nonlinear interaction, E } , so 
that quantum effects are weak (strong) for 



9 



VEJE^ 



(3) 



small (large) compared to 1. The temperature is also 
conveniently measured in the scale of E 3 by the dimen- 
sionless parameter t = T/E 3 . 

For the two-condensate case, starting from the uncer- 
tainty principle (An 2 ) (Asin 2 <^) > (ft 2 /4) (cos<^), the 
authors of Ref. [|3| introduced the coherence factor as the 



thermodynamic average 



a(g,t) = (cos<£) 



(4) 



a parameter that can be directly related with the fringe 
structure in momentum space, which is due to the onset 
of coherence between the weakly linked condensates and 
can be detected by light scattering experiments. Indeed, 
the observed momentum distribution n(p) turns out to 
be described, in terms of the single-condensate distribu- 
tion n (p), as UU nip) = 2 [l+a(g, t) cos(p x d/h)} n (p), 
where d is the distance of the condensate traps (assumed 
along the x axis). In an analogous way, purely quantum- 
mechanical correlations between identical particles lead 
to a reduction of the neutron scattering intensity of pro- 
tons or deuterons which can be explained by a similar 
mechanism |7|. At zero temperature, 1 — a measures the 
decoherence effect of the only zero-point quantum fluctu- 
ations. The value a = 1 occurs when the quantum cou- 
pling 5 = and corresponds to the classical limit, while a 
decreases for increasing values of g. When the tempera- 
ture is finite, a decreases further due to the thermal fluc- 
tuations that contribute to the destruction of coherence. 
This thermal decoherence turned out to be significant 
even at very low temperatures where the condensation 
starts to occur, as shown in Ref. j| for the single BEC 
junction. However, the analogous calculation for an ar- 
ray of condensates, that involves the phase correlation 
function 



G r ig, f ) = ( COS (^i - 0i+r)) , 



(5) 



cannot be performed by a direct numerical solution of 
the Schrodingcr equation, as done in Ref. ||. 

We note that an ideal tool to face such calculations is 
represented by the effective potential method || . Several 
applications H proved indeed its usefulness: in particu- 
lar, it can be used for one- (ID) and two-dimensional 
(2D) arrays JlQ], where, otherwise, the exact quantum 
solution can be obtained only by resorting to much heav- 
ier quantum Monte Carlo simulations. The method is a 



semiclassical expansion in terms of the quantum coupling 
g = \/ ' E c JWj . The only pure-quantum part of the fluc- 
tuations is considered in the self-consistent Gaussian ap- 
proximation so that the quantum Hamiltonian is reduced 
to a classical-like one, where the parameters of the poten- 
tial are renormalized by the pure-quantum effects at any 
temperature. Of course this renormalization vanishes at 
high temperatures, (/3 = \/T — > 0), the approximation 
thus being more and more accurate. 

Let us first consider the two- condensate case. The 
renormalization parameter that takes into account the 
quantum fluctuations only, is related to the pure-quan- 
tum Gaussian spread of the phase; it turns out to be 



<9- 




D o (g,t) = A0l t -A0 2 class 



J 
2 ko 



where f(g,t) = j3huj/2 = gno/(2t) and the renormal 
ized frequency is u) — kqujo; the pure- quantum Hartree 
factor 



- (coth/ ) , (6) FIG. 1. Coherence factor a(g,t) = (cos 0), Eq.(|). Solid 

K0 I line, g = (classical); dotted line, g = 0.5; dashed line, g = 1; 

dash-dotted line, g = \/3. The symbols are the exact result [bj 
for g = 1 (full circles) and g = \/3 (open circles). The latter 
value corresponds to E c /E 1 = 3. 



Ko (g,t) = e- D ^/4 



(7) 



is determined self-consistcntly with Do(g,t), which is a 
decreasing function of temperature and vanishes for t — > 
oo. The effective classical Hamiltonian || corresponding 
to Eq. (|2|) bears the same functional form, 



nn — — — n 



E 3 Ko(g,t) cos <p 



(8) 



with the renormalized Josephson coupling E } K^{g,t) 
(some additive uniform terms, which do not affect the 
calculation of thermal averages, have been neglected). In 
the effective potential formalism any quantum average is 
evaluated by means of a classical-like expression, that for 
the case of the coherence factor reads 



*(#>*) = £- / d<P (ulcosip) exp 
= «o 2b(«o/*) > 



COS(Z> 

t Y 



(9) 



where Xq is the logarithmic derivative of the modified 
Bessel function of the first kind Io(x): 



1 {x) 



dlnJo(x) h(x) 



dx 



Io(x) 



For 5 = one has Ko(0, i) = 1 and Eq. (0) reduces to the 
exact classical coherence factor a c i(t) — Xo(l/t), while 
it is worthwhile to notice that within our approxima- 
tion quantum effects are described by a renormalization 
of the temperature scale. The results are reported for 
different quantum-coupling values in Fig. [j], where com- 
parison with the exact ones is made. It appears that the 
method is reliable as long as the renormalizations are 
small, Do < 1. Note that Do is a decreasing function of 
t, so this condition is verified also at strong coupling if t is 
high enough, while the validity in the whole temperature 
range requires _D (g,0) < 1, i.e., g < 1.5. The result for 



g = 1 is in very good agreement with the corresponding 
one of Ref. [p| in the whole temperature range, while, as 
expected, the result for g = v3 matches the exact data 
for t > 0.5. This testifies to the reliability of the method, 
which can be then pursued also for condensate arrays, 
where quantum fluctuations are expected to be weaker. 

Let us now come to the case of a ID condensate array. 
The effective Hamiltonian corresponding to Eq. (Q) reads 



JV C 



^ = E 



E c 



E 3 K 2 (g,t)cos(ipi - tpi+i) 



(11) 



where now the pure-quantum Hartree factor is given by 
«( ff ,*) = e-^'*)/* ) (12) 

and is to be self-consistently evaluated together with the 
renormalization parameter 



_ . . g 1 ^— v I . k ( . 

Vi(q,t) = — 7= > sin — coth 

Uy ' j 2V2kN„ 4-1 2 V 



t/fc 



1 

fk 



(10) which is a decreasing function of t since 



fk(g,t) 



(3huj k 



gn 

s/2t 



sm ■ 



(13) 



(14) 



indeed, the renormalized dispersion relation takes the 
form hujk = y/2 KU>o\sm(k/2)\. The approximation 
holds if T>\(g,t) < 1, and since at zero temperature 
^i(fl)O) = 9/(V2t{k), the condition amounts to require 
.9 ^ 3.5. 

For the correlation function (|5|) the effective-potential 
recipe gives O] 



G r (g,t) = e 



- p -T> r /2 



'M^/t)Y 



(15) 



1.0 



0.8 



0.6 



<9- 
I 

< J? 0.4 



O 
u 



0.2 



0.0 




„ 6 



0.0 



0.5 



1.0 



1.5 
t 



2.0 



2.5 



3.0 



FIG. 2. Nearest-neighbor correlation function of the ID 
condensate chain Gi(g,t) = (cos(tpi — ipi+i)), Eq.(|17|). Solid 
line, 5 = (classical); dotted line, g = 0.5; short-dashed line, 
g — 1 ; long-dashed line, g — 2; dash-dotted line, g = 3. 



with Zo (#) a s in Eq. (|10|) ; the analytical evaluation uses 
periodic boundary conditions and the translation invari- 
ance of the potential |ll| | . The renormalization parame- 
ters for the rth neighbors read 



zv(ff,*) 



2V2 kN s 



s -^sm 2 (rk/2) 

k 



sin(fc/2) 



coth/ fc - — 

Jk 



(16) 



for r — ► 00 one can replace sin 2 (rfc/2) by its average 1/2 
and get the asymptotic value T> oa {g 1 1), which is logarith- 
mically divergent for t = 0. In particular, one has the 
power-law asymptotic behavior e~' Doo( ' £ ' ,t ^ 2 ~ f , with 
V = 5/[4\/27TK(g, 0)] = T>i{g, 0)/4. Moreover, the zero-i 
correlation function behaves as G r (g,Q) = e~ 



-XV( fl ,0)/2 



r - '' in agreement with the theoretical prediction for the 
low-coupling phase fll2fl . 

For the first-neighbor correlation, we can define the 
analog of the coherence factor, Eqs. (||), which in the 
effective-potential approximation bears the same form of 
Eq. ©, 



(cos(ipi - (f>i+i)) = k 2 l (n 2 /t) ; 



(17) 



it is plotted in Fig. 0. Comparing with the two- 
condensate case (e.g., curves for g = 0.5 and 1), it ap- 
pears that the quantum correction is weaker. It is worth- 
while to notice that in the classical case, i.e., g = 0, 
Eq. ([l5]) reduces to the correct classical expression. 

The formulas above can also be easily generalized [flOfl 
to the 2D case, of recent experimental realization |13f| , 
but the involved classical-like integrals cannot be evalu- 
ated analytically as in Eq. ( |15| ) , so a numerical approach 
is needed, e.g., Monte Carlo (MC) simulation (still much 
easier than quantum MC). 

For the multicondensate array (with spacing d along 
the x direction), the momentum distribution n(p) = 
N s tiq{p) 5(px) exhibits a fringe structure determined by 
the structure function 





p x d/h 

FIG. 3. Structure function E,(p x ;g,t) of the ID condensate 
chain [Eq.jigj)] at fixed t = 0.5. Solid line, g — (classical); 
dotted line, g — 1; dashed line, g — 2; dash-dotted line, 
g — 3. The low- lying curves report the small contribution 
of the term 3i(p x ;g,t) for the same values of g. One single 
periodic interval is displayed. 
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e ~ir Px d/H e -V r /2 



(cos(^> 4 -lfi i+r )) 
[lo(n 2 /t)] 1 



(18) 



where Eq. (|15| ) has been inserted; clearly, 3(p x ) = 
S(p x + 2irh/d) is periodic. In this expression one can 
take advantage of the fact that the renormalization coef- 
ficients T> r (g, t) rapidly converge to the asymptotic value 
T^oo(9jt), separating it as 



E(p x ) = Z (p x ) +Si(p x ) 



-oiPx) 

'-liPx) 



-B»/2 



1 



1 + a 2 — 2a cos(p x d/h) 



E, 



e -ir P *d/h a \r\( e -V r /2 



-Ccx J /2\ 



where the quantity 



*(g,t)=l (K 2 (g,t)/t) 



(19) 
(20) 
(21) 

(22) 



characterizes the coherence effect. Looking at Eq. ( pT| ) 
it appears that the difference in parentheses cuts off the 
contribution from large values of r, so this term does not 
show sharp features and the fringe structure is mostly 
given by Eq. (EQ). In the classical limit (<? = 0) V r — 
vanishes and so does 5i(p x ), and the exact classical 
result H is immediately recovered since a(0,t) = a c \{t). 
Representative results for the full structure function 
E(p x ) and the correction Si(p x ) are reported in Fig. 
for some values of the quantum coupling g at fixed tem- 
perature t = 0.5. The stronger the quantum effects (i.e., 
high g and low t) the higher is Si^), but in the consid- 
ered parameter range the relative contribution of 3i(p x ) 
is small and broad, confirming that the main structure 
is contained in Ho(p x ). In Fig. |J the structure function 
S(p x ) is plotted for fixed g = 0.5 and different temper- 
atures, in a wider range of p x so to display the periodic 
behavior. 
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FIG. 4. Structure function S(pa.; <?,£) of the ID condensate 
chain [Eq.@], at fixed g = 0.5. Solid line, t = 0.2; dotted 
line, t — 0.5; dashed line, t — 1. 



Finally, let us discuss how phase fluctuations do af- 
fect the observable fringe structure. Neglecting the con- 
tribution of Si(p x ), the structure function shows max- 
ima Sm = e~ T '< x l' 2 (1 + a)/(l — a) around the points 
p x = 2£nh/d (£ integer), and minima S m = e~ T '<^l 2 (1 — 
a)/(l + a) in between. The contrast factor is defined as 



Q(5)*) = ( s m - S m )/(H M + Sm) 



(23) 



(Q ~ 1 means good contrast), while the peak width T, 
defined as the full width at half height between maxima 
and minima, is 



T(g, t) = (4h/d) tan" 1 [(1 - a) /{I + a)] 



(24) 



For r <C 2nh/d, i.e., when a is not far from its zero- 
temperature value u{g, 0) = 1, the peaks are sharp 
and have a quasi-Lorentzian shape. However, Eq. ( P0|) 
shows that the overall intensity is depressed by the fac- 
tor e -1300 / 2 , which is plotted in Fig. g together with Q 
and r, for different coupling values: it appears that the 
main effect of quantum fluctuations is to weaken the over- 
all intensity, while the contrast and the peak width are 
predominantly driven by temperature. 
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